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ABSTRACT 

This paper discusses the X-ray halo around the Swift 7-ray burst GRB 050724 (z = 0.258), detected by the 
Swift X-Ray Telescope. The halo, which forms a ring around the fading X-ray source, expands to a radius of 
200" within 8 ks of the burst exactly as expected for small-angle X-ray scattering by Galactic dust along the 
line of sight to a cosmologically distant GRB. The expansion curve and radial profile of the halo constrain the 
scattering dust to be concentrated at a distance of D = 139 ± 9 pc (from Earth) in a cloud/sheet of thickness 
AD < 22 pc. The halo was observed only out to scattering angles of 200", for which the scattering is 
dominated by the largest grains, with a maximum size estimated to be a max w 0.4 — 0.5 /im. The scattering- 
to-extinction ratio was estimated to be T scat /Av J> 0.022; this is a lower limit to the true value because 
contribution from smaller grains, which scatter to larger angles, was not directly observed. The line-of-sight 
to the GRB passes close to the Ophiuchus molecular cloud complex, which provides a plausible site for the 
scattering dust. 

Subject headings: gamma rays: bursts — X-rays: general — Galaxy: structure — ISM:dust 
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1. INTRODUCTION 

Small-angle scattering of X-rays by dust grains can produce 
a 'halo' around a distant X-ray source, with a radial inten- 
sity distribution that depends on the dust properties and loca- 
tion (Overbeck 1965; Hayakawa 1970; Triimper & Schofelder 
1973). This effect was first detected by Rolf (1983) using Ein- 
stein observations of the bright X-ray binary GX 339 — 4, and 
has subsequently been observed around several other bright 
Galactic X-ray sources (e.g. Catura 1983; Predehl & Schmitt 
1995;Predehl et al. 2000). 

Gamma-ray bursts (GRBs) produce high X-ray fluxes for 
short periods, typically <; 1000 s. Viewed through a sub- 
stantial column of Galactic dust, these impulsive X-ray events 
may produce halos that appear to expand on the sky because 
X-rays scattered at larger angles travel a slightly longer path 
length to the observer and suffer an increased time delay. This 
time delay, between the direct and scattered light, means the 
later observations of the scattered light provide a view of 
the GRB X-ray emission at earlier times and can, in princi- 
ple, provide details of the location, spatial distribution and 
properties of the dust (e.g. Alcock & Hatchett 1978; Klose 
1994; Miralda-Escude 1999; Draine & Bond 2004). In par- 
ticular, assuming the GRB to lie at a cosmological distance, 
the distance to the scattering dust can be measured from the 
radial expansion of the halo, because the halo size increases 
as 9 2 — 2tc/D, where r = t — tx is the time since the 
pulse of illuminating X-rays and D is the distance between 
the dust and observer. Observations of GRB dust haloes can 
thereby provide very accurate distances to Galactic structures, 
as first demonstrated by Vaughan et al. (2004) using an XMM- 
Newton observation of the Integral GRB 031203. 

Electronic address: sav2@star.le.ac.uk 

1 X-Ray and Observational Astronomy Group, University of Leicester, Le- 
icester, LE 1 7RH, UK 

2 INAF-Osservatorio Astronomico di Brera, Via Bianchi 46, 23807 Mer- 
ate, Italy. 

3 Department of Astronomy & Astrophysics, Pennsylvania State Univer- 
sity, 525 Davey Lab, University Park, PA 16802, USA. 

4 Universita degli studi di Milano-Bicocca, Dipartimento di Fisica, Piazza 
delle Scienze 3, 1-20126 Milan, Italy. 



This paper discusses the X-ray halo around GRB 050724, 
the first halo discovered by Swift (Romano et al. 2005) and 
only the second ever observed around a GRB. The rest of the 
paper is organised as follows: Section [2] discusses the Swift 
observations and basic data analysis; Section [3]describes the 
X-ray properties of the burst counterpart during the first few 
ks of the observation; Section@]discusses the detailed analysis 
of the X-ray halo images; Section|5]uses the archival ROSAT 
and IRAS data to probe the ISM in the direction of the GRB; 
and finally the implications of the results of this analysis are 
discussed in Section [6] Throughout this paper the quoted er- 
rors correspond to 90% confidence regions unless stated oth- 
erwise. 

2. OBSERVATIONS AND DATA ANALYSIS 

At 2005-07-24 12:34:09 UT the Swift Burst Alert Telescope 
(BAT; Barthelmy 2004) triggered and located GRB 050724 
(Covino et al. 2005). The prompt BAT light curve showed 
a relatively short-lived spike, with a FWHM of w 0.25 s 
(Covino et al. 2005; Krimm et al. 2005), plus an extended 
low-flux tail lasting for at least > 200 s. As discussed by 
Barthelmy et al. (2005), this burst would have been classified 
a short burst if observed by BATSE. The total fluence was 
» 6.3 ± 1.0 x 10~ 7 erg cm -2 (Krimm et al. 2005) over the 
15 — 350 keV band. The spectrum evolves from relatively 
hard during the 1 s peak (photon index of T = 1.71 ± 0.16) 
to much softer (T = 2.5 ± 0.2) at 50 - 150 s after the trigger. 

The spacecraft slewed immediately to the BAT on-board 
calculated position; the Swift X-ray telescope (XRT; Burrows 
et al. 2004, 2005a) began observations at 12:35:22.9 UT (in 
automatic state), only 74 s after the BAT trigger, and de- 
tected a rapidly fading, uncatalogued X-ray source at (J2000) 
RA=16 h 24 m 44.4 s , Dec= -27° 32' 28" (with a 90% confi- 
dence radius of 6"; Barthelmy et al. 2005). This position is 
close to the Galactic plane (I = 350.37°, b = +15.10°); the 
expected column density of neutral Galactic gas along this 
line-of-sight (LoS) is Ah = 1-46 x 10 21 cm -2 (as measured 
from 21 cm maps; Dickey & Lockman 1990), and the redden- 
ing E{B - V) = 0.59 (from the IR dust maps of Schlegel, 
Finkbeiner& Davis 1998). 
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FIG. 1 .— 0.2-5 keV XRT images of GRB 050724 taken in three time intervals covering the first orbit (343 - 2243 s after T ), second orbit (6068 - 8127 s 
after To) and later orbits of Swift, plotted in detector coordinates. The images were smoothed using a Gaussian of width 5 pixels (<r). The X-ray counterpart of 
GRB 050724 is in the centre of each frame. During orbit 1 the halo expands from a radius of 20 to 44 pixels, during orbit 2 it is visible as a faint ring with a 
radius of ~ 85 pixels, and in later orbits the halo is not detected. 



Follow-up observations revealed a variable optical and ra- 
dio source within the XRT error circle (4.3 " from the XRT 
position; Gal- Yam et al. 2005; Soderberg et al. 2005). A 50 ks 
Chandra observation further refined the position of the fading 
X-ray source and showed it to be coincident with the opti- 
cal and radio transient (Burrows et al. 2005b). Prochaska et 
al. (2005) identified the host galaxy of the transient to be a 
massive early-type galaxy at a redshift of z = 0.258 ± 0.002 
(from Ca H+K and G-band absorption). See also Berger et al. 
(2005) for a discussion of the host galaxy optical spectrum. 

The XRT collected data in Windowed-Timing (WT) mode 
from 79 — 342 s post-burst (WT mode allows imaging in one 
spatial dimension only), after which time the data were taken 
in Photon Counting (PC) mode. WT mode allows imaging 
in one spatial dimension only, whereas PC mode produces 
two-dimensional images, see Hill et al. (2004) for a full de- 
scription of the XRT operating modes. The XRT data were 
processed by the Swift Data Center at NASA/Goddard Space 
Flight Center (GSFC) to level 1 data products (calibrated, 
quality flagged event lists). These were further processed with 
xrtpipeline vO . 8 . 8 into level 2 data products. In the 
subsequent analysis only event grades — 2 were used for the 
WT mode data and only grades — 12 were used for the PC 
mode data. 

3. EARLY X-RAY EMISSION FROM GRB 050724 

Spectra and light curves were extracted from the WT data 
using a 40 pixel (94") wide box to define the source and back- 
ground regions. For the PC mode data, source counts were ac- 
cumulated from within a circle of radius 25 pixels (59"), and 
background data were accumulated within an annulus having 
inner and outer radii of 60 and 100 pixels (141 — 236"). The 
appropriate ancillary response matrices were generated using 
xrtmkarf vO.4.14. 

The spectrum of the WT data (79 — 342 s post-burst) was 
binned such that each spectral bin contained at least 20 counts 
and fitted over the 0.6 - 10 keV range 5 using XSPEC (Arnaud 
1996). The spectral model comprised a power law contin- 
uum modified by two neutral absorbers (modelled using the 

5 The bandpass for the WT-mode spectral fit was restricted to energies 
above 0.6 keV to avoid calibration uncertainties at lower energies. 
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FIG. 2. — Ax 2 contours for absorption in the host Galaxy (z = 0.258) 
against Galactic absorption (z = 0.0). The contours correspond to A\ 2 = 
2.30, 4.61 and 9.21, respectively. The dotted line indicates the expected /Vh 
based on the 21 cm measurements of Dickey & Lockman (1990). 



TBabs code of Wilms, Allen & McCray 2000), one at z = 
0.0 and the other at z = 0.258, corresponding to neutral inter- 
stellar gas in our Galaxy and intrinsic to the GRB host galaxy, 
respectively. Allowing both column densities to be free pa- 
rameters, this simple model gave a very good fit (y 2 = 355.89 
for 334 degrees of freedom, dof) with best-fitting parameters 
as follows: T = 1.94 ± 0.05, N H {z = 0.0) = 5.9±?j x 10 21 
cm" 2 and iV H (z = 0.258) < 2.4 x 10 21 cm" 2 . The Galactic 
column density inferred from the X-ray spectrum exceeds that 
expected based on the 21 cm measurements (see Section |3, 
by 4.4^'g x 10 21 cm~ 2 , and yet there is no clear excess ab- 
sorption intrinsic to the host galaxy. Figure [2] shows the Ay 2 
contours for the two absorption column densities. 

4. X-RAY HALO ANALYSIS 

4.1. Expansion curve 

Figure ^ shows the XRT (PC mode) images from the first, 
second and later orbits of Swift, over the 0.2 — 5 keV en- 
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FIG. 3. — Radial profile of counts about the GRB 050724 afterglow from 
five time intervals. The times are marked in each frame. The arrows mark the 
maximum likelihood estimate of the halo radius, 8, derived from fitting these 
profiles, and the horizontal bar indicates the 68.3% (lcr) confidence interval 
on the radius. 
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FIG. 4. — Change in C-statistic as a function of halo position. The min- 
imum of each curve represents the maximum likelihood estimate of the halo 
radius. Nominal 68.3 and 90% confidence regions can be defined by the 
intervals bounded by AC = 1.0 and 2.706, respectively. 



ergy range (only events with energies 0.2 — 5 keV were 
use d sin ce the halo has a soft spectrum, as discussed in Sec- 
tion 14.31 and was not detected above 5 keV). The first orbit 
showed extended emission around the point-like GRB which 
was dispersed in the second orbit and was difficult to detect 
in later orbits. In order to better quantify this extension, im- 
ages were accumulated in five non-overlapping time intervals: 
343 - 543, 543 - 943, 943 - 1543 and 1543 - 2243 s from 
the first orbit and 6068 — 8127 s from the second orbit. 

Radial profiles were calculated for each image by accumu- 
lating the counts in annuli centred on the GRB (see Figure[3}. 
The expanding halo was revealed by the small but significant 
localised excess of counts, beyond the unresolved GRB 6 . The 
radius of the halo as a function of time, 0(t), was estimated by 
fitting each radial profile (using XSPEC) with a model com- 
prising the PSF of the source (a King profile with parameters 
taken from ground based calibration data; Moretti et al. 2004, 
2005) plus a constant background level (per pixel), integrated 
over the annular bin. A Gaussian was added to represent the 
profile of the halo, fixing the width to be equal to 10% of the 
radius: ag = 0.16 (see Section l4~2l >. The C-statistic was used 
to find the maximum likelihood model parameters, appropri- 
ate for the case of few counts per bin (C — — 2 log C, where C 
is the likelihood; Cash 1979). For each profile, the C-statistic 
was minimised and the maximum likelihood position of the 
Gaussian was used as the radius of the halo in that image. The 
90% confidence region for the radius was estimated using the 
AC values (see Figure^). 

Figure |3 shows the radius of the halo in each of the five 
time bins, as measured using the radial profiles. These data 
were fitted with the function 6{t) = N(t - tx) 1/2 , which 
describes the radial expansion of a dust-scattered halo fol- 
lowing a (short) X-ray pulse at time tx (e.g. Triimper & 
Schofelder 1973; Klose 1994; Miralda-Escude 1999; Draine 
& Bond 2004). Strictly speaking, this assumes a very nar- 
row pulse of X-rays and an intrinsically narrow dust scatter- 
ing screen. The value of tx corresponds to the typical arrival 

6 The Point Spread Function (PSF) of the XRT has a half-energy radius of 
8.8", or FWHMsi 9.6" HEm 1.5 keV (Moretti et al. 2005; Burrows et 
al. 2004) 



time of direct (un-scattered) photons, which will be after To 
since the GRB 050724 remained X-ray bright from To until 
« To 4- 300 s, meaning almost all the direct X-rays actually 
arrived after To. The best-fitting values of the free parame- 
ters were: N = 2.46 ± 0.08 (with 6 measured in arcsec) and 
tx = 145l 6 4 s after T . This model gave a good fit, with 
X 2 = 6.70 for 4 dof. In order to check whether the best- 
fitting value of tx is plausible, the mean arrival time for the 
direct X-rays was estimated using the X-ray light curve from 
the joint BAT-XRT data (Figure 3 of Barthelmy et al. 2005): 
the flux-weighted photon arrival time 7 over this light curve 
was To + 86 s, consistent with the best-fitting tx from the 
radial expansion curve. 

The normalisation, 7Y, describes how fast the halo expands 
and is uniquely determined by the distance to the dust 8 . As- 
suming the GRB to be at a cosmological distance N = 
(2c/ D) 1 / 2 . The best-fitting normalisation predicts a distance 
to the scattering medium of D = 139 ± 9 pc, where the 90% 
confidence region was calculated with tx as a free parameter 9 

This distance measurement is quite robust to the choice of 
function used to parameterise the halo radial profile. This was 
demonstrated by repeating the radial profile fitting using dif- 
ferent functions to model the halo profile. For each one the 
halo expansion curve was measured and the distance calcu- 
lated from its best-fitting normalisation. A Gaussian of fixed 
width (er = 5") gave D = 145 ± 8 pc, a Lorentzian of fixed 
width (r = 10") gave D = 142 ± 8 pc and a King model 
as used to model the PSF (see above) gave D = 144 ± 6 pc. 
These different estimates are all well within the confidence re- 
gion of the original estimate, demonstrating that any system- 
atic errors introduced by the choice of radial profile model are 

7 The flux- weighted photon arrival time was calculated as (tx) = 
J /x(*)t dt/ J fx(t) dt where /x(*) is the X-ray light curve. 

8 The halo shape is affected by the geometry of the scattering region, and 
it may not be circular and centred on the GRB image if the scatterer is in 
the form of a plane-parallel slab not perpendicular to the line-of-sight, or is 
curved (see Tylenda 2004 for a discussion). However, because of the large 
ratio of dust distance to halo size, only the most pathological geometries (e.g. 
slab included at > 89°) will produce any observable effect. 

9 This distance estimate is smaller than that given in Romano et al. (2005) 
because that earlier analysis assumed tx = To , an assumption not made in 
the present analysis. 
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FIG. 5. — Expansion of the ring around GRB 050724 with time. The radii 
were measured by fitting the local maximum in the radial profiles (Fig.l3l. 
The expansion was fitted with the function 9{t) = N(t — fx) 1 ' 2 (marked 
with the dotted line). 
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FIG. 6. — Co-moving radial profile around GRB 050724. This shows the 
counts accumulated in annuli around the GRB summed over the five time 
intervals of Figure|5] after correcting for the expansion of the halo (matching 
the radii to those from the first image). The halo occurs at 8 ?s 45 arcsec 
from the GRB (9 = 0) and is a relatively narrow (ag /8 fa 0.07) structure. 



smaller than the quoted statistical error. 

4.2. Co-moving radial profile 

The radial profile of the halo was examined in more detail 
by summing the five radial profiles (Figure [3j, after adjust- 
ing the size of the radial bins to correct for the expansion of 
the halo. A constant background level, estimated from the 
mean level at large radii, was subtracted from each profile be- 
fore combining them. The result is shown in Figure [5] This 
shows that the average radial profile of the halo, after correct- 
ing for its expansion, was reasonably narrow. Fitting this pro- 
file using a Gaussian to model the halo gave a relative width 
of FWHM(0)/0 ps 0.16 - 0.23 (90% confidence limit). 

The observed radial width of the halo is a combination of 
four terms: (i) the time profile of the illuminating X-ray pulse; 
(ii) the radial distribution of the dust along the line of sight, 
(Hi) the PSF of the XRT, and (iv) the radial expansion of the 
halo over the length of an exposure. Term (i) is negligible 
since the time duration of the initial pulse is short compared 
to the times (since burst) of the halo images (343 — 8127 s). 
The XRT PSF (with FWHM of 9.6") gives rise to a fractional 
width of 20% in the first image to 5% in the last, meaning term 
(Hi) must contribute significantly to the measured width. The 
co-moving radial profile is dominated by counts from later 
images (with longer exposures), during which the fractional 
increase in the halo radius - term (iv) - was ~ 15%. The 
measured width of the halo may therefore be explained en- 
tirely by terms (Hi) and (iv), implying the remaining term, 
the radial thickness of the dust, must be small. A very conser- 
vative limit on the physical thickness of the dust cloud along 
this line of sight is AD ps D x FWHM(0) /0 < 22 pc. Most 
likely the dust is confined in a much narrower region. 

4.3. Halo spectrum 

The spectrum of the halo was extracted from the first or- 
bit of data using an annular extraction region with inner and 
outer radii of 20 and 45 pixels (47 — 106 arcsec), respectively. 
A background spectrum was extracted from an annulus with 
radii of 60 and 100 pixels and an appropriate ancillary re- 
sponse file was generated. Due to the small number of counts 
(ps 100) the spectrum was fitted using the C-statistic, rather 
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FIG. 7. — XRT spectrum of the halo during the first Swift orbit. Crosses 
mark the data and the histogram marks the best-fitting absorbed power law 
model. The data were binned for illustration purposes only. 



than binning to 20 counts per bin and fitting using the x 2 
statistic, which would result in very few bins. The spectral 
model used was an absorbed power law, with the neutral ab- 
sorption column density fixed at the value obtained from the 
early XRT spectrum (iVn = 5.9 x 10 21 cm~ 2 ; section^. The 
data are shown in Figure^ The best fitting photon index was 
r = 2.7lg 3, significantly steeper than the GRB X-ray spec- 
trum (T = 1.94 ± 0.05), as expected due to the strong energy 
dependence of the dust scattering cross section. The model 
gave a good fit, with a rejection probability of only p = 0.31 
estimated using 10 3 Monte Carlo simulations. The mean pho- 
ton energy in the binned, background subtracted halo spec- 
trum was (E) = 1.54 keV. 

4.4. Halo flux as a function of angle 

The flux decay of the halo was also estimated using the ra- 
dial profile of the combined image from orbits 1 and 2 (over 
which the halo is well-detected). Counts were accumulated 
in annuli of width A0 = 17.6" (7.5 pixels) and background 
subtracted assuming a constant background level per pixel es- 
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timated from an annulus enclosing 280 — 370". The sampling 
of the halo as a function of angle is incomplete, due to the 
incomplete time coverage of the observation, therefore radii 
over which the halo was not observed were masked out of the 
radial profile. The masking was done by converting the times 
of the observations into radii using the best-fitting expansion 
curve (Figure and removing all radial bins containing un- 
observed angles. The result is shown in Figure|H] 

The total fluence of scattered light (in the halo), Fh, should 
be related to the direct light (i.e. un-scattered light from the 
GRB), F x , by: 



F H 



Fa + Fx 



= 1 



(1) 



where r scat = aN s is the scattering optical depth, with a 
the total scattering cross-section and jV g is the dust grain col- 
umn density along the LoS (e.g. equation 14 of Mauche & 
Gorenstein 1986, henceforth MG86). For low optical depths 
(7"scat -C 1), where multiple scattering is insignificant, this 
can be approximated by Fh w T scat Fx. The above quantities 
represent the entire halo (integrated over all scattering angles) 
but the scattering cross-section is a function of angle which 
defines the halo fluence per unit angle. In the Rayleigh-Gans 
approximation 10 the differential scattering cross-section for a 
spherical grain of size a (in units of /im) is given by: 



da 



A K a 6 



(1 + cos 2 9) 



(2) 



(equation 2.2 of Hayakawa 1970) where Ae is a normalisa- 
tion that depends on the energy of X-rays being scattered, 
x = (47ra/A) sin(0/2) and ji{x) — (sinx)/a; 2 — (cosx)/x 
is the spherical Bessel function of the first order. The cen- 
tral core of this function is approximately Gaussian, leading 



to a halo with a characteristic size 6 TI 



62A(aEY 



arcsec 



(equation 6 of MG86; Hayakawa 1970). Using equation|2]the 
scattered fluence in the halo per unit scattering angle is: 



dF H 

de 



-A E a 6 2Tr6 



JiO) 



(1 + cos 2 6)N S F X 



<xa b 9 



1 2 



(1 + cos 2 



(3) 



which peaks at 9 » 9 lms . 

This function was fitted to the data (solid curve in Figure|8j 
yielding a best-fitting value of a = 0.367 /im (with a 90% 
confidence region of 0.316 — 0.430 /im), assuming a typical 
ene rgy o f E « 1.54 keV for the scattered X-rays (see sec- 
tion The scattering of soft X-rays is dominated by the 
largest grains, due to the strong dependence of the scattering 
cross-section on grain size (the total scattering cross section, 
integrated over all angles is a oc a 4 ), and so this estimate for a 
corresponds to grains close to the maximum grain size. This 
best-fitting model was integrated over all angles to estimate 
the total halo fluence of Fh ~ 342 ± 34 counts (where the 
la error was approximated by the 10% fractional error on the 
total observed halo fluence). 

The above calculation is valid for a single grain size which 
should be a reasonable approximation given the strong de- 
pendence of the scattering cross-section on a. The presence 

10 Smith & Dwek (1998) demonstrated that the Rayleigh-Gans theory and 
the more exact Mie theory are in close agreement for photons energies E ss 
2 keV and typical grain sizes. 
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FIG. 8. — Halo flux as a function of scattering angle. The crosses show the 
background-subtracted counts per (equal size) radial bin over which the halo 
was detected. The solid curve marks the best-fitting theoretical scattering 
curve dFn/d9 (Equation|3} assuming a single grain size (a = 0.367 /im). 
The dot-dashed line shows the best-fitting function assuming a MRN distri- 
bution of grain sizes (with a max = 0.524 /im). 



of smaller grains will however cause the halo fluence to fall 
off less rapidly at large angles and this extended 'tail' could 
contain substantial additional fluence. In order to investigate 
the consequence of a continuous grain size distribution, equa- 
tion[3]was re-calculated at each angle by integrating over the 
grain size distribution N(a). The size distribution was as- 
sumed to follow that of Mathis, Rumpl & Nordsieck (1977; 
hereafter MRN): 

f N a~ q a min < a < a max 
N(a) = i (4) 
1 elsewhere 

The revised function was fitted to the data assuming a m [ n » 
0.04 /im and q = 3.5, yielding a best-fitting a max = 
0.524 pm (with a 90% confidence region 0.426 - 0.790 /im), 
this is slightly larger than, but not inconsistent with other esti- 
mates derived from X-ray haloes (MG86; Predehl et al. 1991; 
Clark 2004). The best-fitting model is also shown in Figure[S] 
(dot-dashed curve). The fit is not very sensitive to flmin as 
long as a max /a m ; n ;> 5 and q J> 3.5. The total fluence in 
this model, integrated over all angles, was virtually identical 
to the single grain size model (w 343 ± 34 counts), although 
the fluence could be higher if there is a greater contribution 
from smaller grains scattering to larger angles that were not 
directly observed. 



4.5. Scattering optical depth 

Given the above estimate of the halo fluence, Fh, it is pos- 
sible to estimate the dust scattering optical depth once the 
direct GRB X-ray fluence, Fx, is known. This latter quan- 
tity can be estimated by integrating the combined XRT-BAT 
light curve of the burst. Figure 3 of Barthelmy et al. (2005) 
shows the X-ray light curve of GRB 050724 from T . This 
was produced using the XRT observations after Tq + 79 s, 
and the 15 — 25 keV BAT light curve extrapolated into the 
XRT band based on a spectral fit to the simultaneous BAT and 
XRT spectrum. During the first ~ 300 s after the trigger the 

" Using a m j n = 0.004 /im or 0.1 fim changed the fit statistic by only 
Ax 2 = +0.009 or +0.038, respectively. 
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X-ray source was bright; the WT data contain ~ 1.18 x 10 4 
source counts (0.2 — 10 keV). The extrapolated BAT data pre- 
dict almost exactly the same number of counts prior to the 
start of the XRT/WT observation. Although there will be a 
systematic error associated with the extrapolation of the BAT 
data to lower energies, the tight agreement between BAT and 
XRT data in the interval of overlap suggests this is likely to 
be small (Barthelmy et al. 2005). Together these give a total 
of « 2.42 x 1 4 counts. Combining this with the Fn value 
of section l4~4l gave r scat « Fn/Fx ~ 0.014. Making an ap- 
proximate correction to an effective scattering optical depth at 
1 keV gives r scat (l keV) m T scat (E)E 2 = 0.034 (with E the 
mean scattered photon energy in keV). 

5. ISM TOWARDS GRB 050724 

Archival ROSAT and IRAS data were used to constrain the 
ISM in the direction of GRB 050724. The top two panels 
of Figure [9] show the soft X-ray sky around GRB 050724, 
derived from the ROSAT All-Sky Survey (RASS) data in the 
0.4 — 1.2 keV band, and the bottom panel shows the IRAS All- 
Sky Survey 100 /im map of the region. These images reveal 
structures dominated by the Ophiuchus molecular cloud com- 
plex casting the dark 'shadows' in the ROSAT map (middle 
panel) and enhanced IR dust emission in the IRAS map (lower 
panel). The main pOph molecular cloud (L1688;Klose 1986) 
is the in the bright region in the upper-left part of the IRAS 
map. These data were used to estimate the column density in 
the cloud by comparing the column density estimates for the 
GRB position to those at locations A' and 'B,' chosen to be 
representative of nearby but off-cloud lines-of-sight. 

The column densities were estimated in three independent 
ways, and the results are presented in Tabled The first esti- 
mates comes from the Dickey & Lockman (1990) 21 cm maps 
of TYji. This shows little change in N-r between the three posi- 
tions (see Tabled, and this was confirmed by the 21 cm map 
of the region shown by de Geus & Burton (1991; see their fig- 
ure 5). The column density due to the intervening cloud in the 
direction of GRB 050724 was estimated from the difference 
between the column density towards the GRB and the mean of 
that measured for A and B. This gave N H = 0.05±0.20xl0 21 
cm~ 2 (where the uncertainty is a very rough estimate based 
on half the difference between A and B). In other words there 
is no excess H I in the direction of GRB 050724 due to the 
intervening dust cloud. 

The second estimate comes from the IR reddening maps 
from Schlegel et al. (1998) This gave E(B - V) values, 
based on the COBE temperature and IRAS intensity at the 
position of interest, of 0.59, 0.17 and 0.25 for the GRB, 
A and B, respectively. These correspond to optical estinc- 
tions of Ay = 2.37, 0.68 and 1.01, respectively, assuming 
R v = A V /E(B - V) = 4.0 as measured by Vrba, Coyne 
& Tapia (1993) in the direction of the p Oph cloud. The red- 
dening values were converted to hydrogen column densities 
assuming 7V H = 4.93 x 10 21 E(B - V) cm" 2 (Diplas & Sav- 
age 1994). In this case the excess column density over the 
GRB was jV H = 1-9 ± 0.2 x 10 21 cm -2 . However, as noted 
by Jura (1980), Vuong et al. (2003) and others, the dust-to-gas 
conversion factors in the direction of the Oph cloud are uncer- 
tain due to the complexity of modelling the cloud medium. 

The third estimate was derived from RASS spectra ex- 
tracted 12 using 0.5° radius regions centred on the positions 

12 ROSAT count rates in different channels were extracted using the 
HEASARC "X-ray background tool" available via the following address: 




FIG. 9. — X-ray and IR images of the sky around GRB 050724. Upper 
panel: 0.4- 1.2 keV ROSAT All-Sky Survey (RASS) map in Galactic coordi- 
nates spanning 102° on a side. The position of GRB 050724 is marked with 
the small filled circle, the large circle marks the extent of the super-bubble 
centred on the Sco-Cen OB association. Middle panel: 10° X 10° close up 
of the RASS map around GRB 050724 (marked with the filled circle). The 
dotted circle has a radius of 0.5°, which corresponds to fs 2.4 pc at the dis- 
tance to the scattering dust. Lower panel: 100 fim IRAS map on the same 
scale as the middle panel. Clearly the IR emission and soft X-ray absorption 
are well correlated. Again the position of the GRB is indicated by the filled 
circle. Positions A' and 'B ' mark the off-cloud regions discussed in the text. 
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TABLE 1 

Compilation of Galactic column density estimates 



Technique Reference N H [GRB] 7V H [A] 


N H [B] 




(10 21 cm" 2 ) (10 21 cm" 


" 2 ) (10 21 cm" 2 ) 


21 cm 1 


1.46 1.22 


1.61 


IR map 2 


2.92 0.83 


1.24 


RASS 3 


3.57 0.52 


1.28 



REFERENCES: (1) Dickey & Lockman (1990), (2) Schlegel et al. (1998), 
(3) this work 



of the GRB, A and B. These were fitted with the model de- 
veloped by Willingale et al. (2003) but, due to the complex- 
ity of the model and low resolution of the data, most pa- 
rameters were fixed at the values found from an analysis of 
nearby XMM-Newton fields (M. A. Supper, in prep.). The 
two parameters allowed to vary in the fit were the foreground 
absorption column density and the emission measure of the 
"Loop I super-bubble" (indicated by the large circle in the 
top panel of Figure [9jl. The excess column density towards 
GRB 050724 was estimated, by comparison with A and B, 
at A H = 2.7 ± 0.4 x 10 21 cm~ 2 . Given the systematic er- 
rors involved in modelling the complex X-ray spectrum, and 
converting from dust reddening factor to equivalent hydrogen 
column density, the X-ray and IR derived estimates are in rea- 
sonable agreement. 

6. DISCUSSION 

The X-ray halo around GRB 050724, as observed by the 
Swift XRT, has provided accurate information on the Galac- 
tic dust distribution in this direction. In particular, the narrow 
halo must be caused by a concentration of dust at a distance of 
D = 139 ± 9 pc from Earth, in a cloud or sheet with a thick- 
ness AD < 22 pc. The radial profile of the halo constrains 
the size of the largest grains to be a max ~ 0.4 — 0.5 pm, 
but the contribution from smaller grains is unconstrained be- 
cause the halo was not detected at scattering angles larger 
than 9 w 200". The estimated value of a max is higher than 
the often used value of 0.25 pm (MRN; see also Predahl & 
Schmitt 1995), although this is not without precedent as Jura 
(1980) previously noted that the non-standard optical extinc- 
tion curves in the direction of p Ophiuchus may signify the 
presence of unusually large grains. 

The LoS to GRB 050724 includes the Upper Scorpius sub- 
group of the Scorpius-Centaurus OB association at a mean 
distance of 145 ± 2 pc (derived from HIPPARCOS data; de 
Zeeuw et al. 1999). There is considerable structure in the 
extinction maps around this region, including material be- 
longing to the Ophiuchus molecular cloud complex (de Geus 
1992). Knude & H0g (1998) detected a sharp rise in the red- 
dening at 120 — 150 pc in the direction of the p Oph star form- 
ing region, which is consistent with the distance measured for 
the dust responsible for the X-ray halo. If the X-ray scatter- 
ing dust is associated with part of the Oph molecular cloud 
complex, this gives the most accurate measurements to date 
for the distance to the cloud and its physical thickness. 

The excess optical extinction in the direction of the 
GRB 050724, compared to the mean of the off-cloud positions 

|http : / /heasarc ■ qsf c ■ nasa ■ qov/| 



A' and 'B', was Ay — 1.5 mag (Section^. In combin ation 
with the estimated halo scattering optical depth (Section ^. 5> 
this gives T scat /Av ~ 0.022 at 1 keV. This is lower than the 
value of 0.056 estimated by Predehl & Schmitt (1995) using 
ROSAT observations of bright Galactic X-ray sources, and a 
factor of a few lower than the value derived from the model of 
Draine (2003). However, the halo around GRB 050724 was 
only observed at small angles (8 <J 200") and so the flux con- 
tribution from smaller grains is largely unconstrained, mean- 
ing that the true halo flux, and hence T scat , is probably larger. 

As shown by Figure [9] the X-ray absorption and IR dust 
emission are well-correlated, demonstrating they are caused 
by the same medium. By contrast the 21 cm map shows little 
correlation with these images (Figure 5 of de Geus & Bur- 
ton 1991), suggesting the cloud medium has a comparatively 
low density of atomic H I. The most obvious explanation is 
that the hydrogen is molecular, although there is no obvious 
CO emission from this location in figure 3 of de Geus, Bronf- 
man & Thaddeus (1990). The Ly-a map of the region (see 
figure 13 of de Geus et al. 1990) shows two nearby H II re- 
gions, known as S9 and RCW 129, around the stars a Sco 
(Bl III) and r Sco (B0 V), respectively. It is therefore also 
plausible that some fraction of the atomic hydrogen is ionised 
by these nearby young stars (with distances of ~ 140 pc; 
Shull & van Steenberg 1985). The total Galactic hydrogen 
column density along the LoS to GRB 050724 is therefore 
the sum of that revealed by the 21 cm H I maps and that 
inferred from the molecular cloud (mostly H2 and/or H II). 
The 21 cm measurement of the total H I column density gave 
1.5 x 10 21 cm~ 2 , whereas the column density in the dust 
cloud was measured at 1.9 ± 0.2 x 10 21 cm -2 from the ex- 
cess IR emission and 2.7 ± 0.4 x 10 21 cm -2 from the excess 
soft X-ray absorption (Section [5}- The total hydrogen col- 
umn density is therefore in the range 3.4 — 4.2 x 10 21 cm -2 , 
slightly lower than that measured from the GRB X-ray spec- 
trum (4.4 — 6.2 x 10 21 cm~ 2 ; Section QJ. It is of course 
possible that this extra X-ray absorption along the LoS to 
GRB 050724 (AA H ~ 1.5 x 10 21 cm~ 2 ) is caused by either a 
modest column of molecular or ionised gas on the far side of 
the Sco-Cen super-bubble, or due to cold gas in the GRB host 
galaxy, the formal constraint on which was Ah < 2.4 x 10 21 
cm~ 2 (Section[3}. 

Note added in proof: Since the completion of this paper 
Tiengo & Mereghetti (2005) have announced the detection of 
another dust-scattered X-ray halo around a GRB. This brings 
the total number of known GRB X-ray halos to three. 
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